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ABSTRACT 

We analyze the dust attenuation properties of a volume-limited, optically- 
selected sample of normal star forming galaxies in nearby clusters as observed 
by GALEX. The internal attenuation is estimated using three independent in- 
dicators, namely: the ratio of the total infrared to far-ultraviolet emission, the 
ultraviolet spectral slope jS and the Balmer decrement. We confirm that nor- 
mal galaxies follow a Ltir/ Lpuv — P relation offset from the one observed for 
starburst galaxies. This offset is found to weakly correlate with the birthrate 
parameter, thus with the galaxy star formation history. We study the correla- 
tions of dust attenuation with other global properties, such as the metallicity, 
dynamical mass, ionized gas attenuation. Ha emission and mass surface den- 
sity. Metal-rich, massive galaxies expected, more heavily extinguished in 
the UV than are small systems. For the same gas metallicity normal galaxies 
have lower Ltir/ Lpjjy ratio than starbursts, in agreement with the difference 
observed in the Ltir/ Lpuv — P relation. Unexpectedly we find however that 
normal star forming galaxies follow exactly the same relationship between metal- 
licity and ultraviolet spectral slope (3 determined for starbursts, complicating 
our understanding of dust properties. This result might indicate a different dust 
geometry between normal galaxies and starbursts, but it could also be due to 
aperture effects eventually present in the lUE starbursts dataset. 
The present multiwavelength study allows us to provide some empirical rela- 
tions from which the total infrared to far ultraviolet ratio [Ltir/ Lpjjy) can be 
estimated when far infrared data are absent. 
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burst, spirals 

1. Introduction 

The presence of dust in galaxies represents one of the major obstacles complicating 
a direct quantification of the star formation activity in local and high redshift galaxies. 
Absorption by dust grains reddens the spectra at short wavelengths and modifies altogether 
the spectral energy distribution of galaxies. Since the UV radiation is emitted by young stars 
(t < 10^ yr) that are generally more affected by attenuation from surrounding dust clouds 
than older stellar populations, rest-frame UV observations can lead to incomplete and/or 
biased reconstructions of the star formation activity and star formation history of galaxies 
affected by dust absorption, unless proper corrections are applied. 

In recent years our understanding of dust attenuation received a tremendous impulse from 
studies of local starburst galaxies (i.e.Calzetti et al. 1994; Heckman et al. 1998; Meurer et al. 
1999; Calzetti 2001; Chariot & Fall 2000), that were based on three indicators: the ratio 
of the total infrared to far- ultraviolet emission {Ltir/ Lpuv), the ultraviolet spectral slope 
(3 (determined from a power-law fit of the form / ~ to the UV continuum spectrum in 
the range 1300 and 2600 A, Calzetti et al. 1994) and the Balmer decrement. The total- 
IR (TIR) to UV luminosity ratio method (i.e. Buat 1992; Xu & Buat 1995; Meurer et al. 
1995, 1999) is based on the assumption that a fraction of photons emitted by stars and 
gas are absorbed by the dust. The dust heats up and subsequently re-emits the energy 
in the mid- and far-infrared. The amount of UV attenuation can thus be quantified by 
means of an energy balance. This method is considered the most rehable estimator of the 
dust attenuation in star-forming galaxies because it is almost completely independent of the 
assumed extinction mechanisms (i.e. dust/star geometry, extinction law, see Buat & Xu 
1996; Meurer et al. 1999; Gordon et al. 2000; Witt & Gordon 2000). When the spectrum is 
dominated by a young stellar population the ultraviolet spectral slope is found to have 
a weak dependence on metallicity, IMF, and star formation history (Leitherer & Heckman 
1995). Thus the difference between the observed (3 and the one predicted by models can 
be entirely ascribed to dust attenuation (Meurer et al. 1999). However in systems with no 
or mild star formation activity the UV spectral slope can be strongly contaminated by the 
old stellar populations, whose contribution increases (3 (flattens the UV continuum, Boissier 
et al. 2005). Thus the spectral slope of mildly star forming systems could be intrinsically 
different from the one of starburst galaxies, even in the absence of dust attenuation (Kong 
et al. 2004). 

Meurer et al. (1999) have shown that in starburst galaxies the total far-infrared to ultraviolet 
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luminosity ratio correlates with the ultraviolet spectral slope, (5 (commonly referred to as 
the IRX-UV relation). They pointed out that this relation allows reliable estimates of the 
attenuation by dust at ultraviolet wavelengths based on 13. 

The Balmer decrement gives an estimate of the attenuation of ionized gas and not of the 
stellar continuum as in the previous two methods. It is based on the comparison of the 
observed Hq;/H/5 ratio with its predicted value (2.86 for case B recombination, assuming 
an electronic density rig < 10^ cm~^ and temperature ~ 10^ K; e.g., Osterbrock 1989). 
Calzetti et al. (1994) found a significant correlation between the ultraviolet spectral slope 
(5 and the Balmer decrement Hq;/H/3. Starting from this empirical relation they obtained 
an attenuation law (known as the Calzetti attenuation law) often adopted to correct UV 
observations for dust attenuation in absence of both far-infrared observations and estimates 
of the ultraviolet spectral slope (Steidel et al. 1999; Glazcbrook et al. 1999). 
Unfortunately the above empirical relations have been established only for starburst galaxies 
and they seem not to hold for normal star forming galaxies. Recently, Bell (2002) suggested 
that quiescent galaxies deviate from the IRX-UV relation of starburst galaxies, because they 
tend to have redder ultraviolet spectra at fixed total far-infrared to ultraviolet luminosity 
ratio. Kong et al. (2004) confirmed this result and interpreted the different behaviour of 
starbursts and normal galaxies as due to a difference in the star formation histories. They 
proposed that the offset from the starburst IRX-UV relation can be predicted using the 
birthrate parameter h (e.g. the ratio of the current to the mean past star formation activity). 
However an independent observational confirmation of the correlation between the distance 
from the starburst IRX-UV relation and the birthrate parameter has not been obtained so far 
(Seibert et al. 2005). Even the Calzetti law does not seem to be universal. Buat et al. (2002) 
showed that for normal star forming galaxies the attenuation derived from the Calzetti law 
is ~0.6 mag larger than the one computed from the Ffir/Fuv ratio and their result has been 
recently confirmed by Laird et al. (2005). 

Why do normal star-forming galaxies behave differently from starbursts? Do normal galaxies 
follow different empirical relations that can be exploited to correct for dust attenuation in 
absence of far infrared observations? If this is the case, is there a transition between starburst 
and normal galaxies? Which physical parameters drive it? Answering these questions will be 
important for a better understanding of the interaction of dust and radiation specifically in 
nearby dusty star forming galaxies, but it also has direct consequences for our understanding 
and interpretation of galaxy evolution in a general context. Firstly it seems mandatory to 
characterize the dust attenuation properties of normal galaxies, to compare them with the 
ones of starbursts and to derive new recipes for the UV dust attenuation correction. This 
topic came once again to the fore with the launch of the Galaxy Evolution Explorer (GALEX). 
This satellite is delivering to the community an unprecedented amount of UV data on local 
and high redshift galaxies that require corrections for dust attenuation but currently lack 
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far-infrared rest-frame data. The time is ripe to explore new methods for correction of 
these data, that might provide new insights on galaxy evolution. Whenever they can be 
combined with other data, GALEX observations provide the best available ultraviolet data 
for studying the dust attenuation properties of galaxies. Multiwavelength photometric and 
spectroscopic observations are in fact mandatory in order to: determine metallicity, ionized 
gas attenuation {A{Ha)), luminosity and mass, test the validity of the relations followed 
by starbursts (Heckman et al. 1998), explore relations that might prove useful to correct 
ultraviolet magnitudes and to compare them with various models of dust attenuation. Recent 
extensive spectroscopic and photometric surveys, hke the Sloan Digital Sky Survey (SDSS, 
Abazajian et al. 2005) and the Two Degree Field Galaxy Redshift Survey (2dF, CoUess et al. 
2001) have opened the path to studies of fundamental physical parameters based on enormous 
datascts. However, spectroscopic observations of nearby galaxies suffer from strong aperture 
effects, making these datasets not ideal for the purpose of the present investigation. In fact, 
Kewley et al. (2005) have recently shown that aperture effects produce both systematic and 
random errors on the estimate of star-formation, metallicity and attenuation. To reduce at 
least the systematic effects they suggest selecting only samples with fibres that capture > 
20% of the hght. This requires z > 0.04 and z > 0.06 for SDSS and 2dF respectively: too 
distant to detect both giant and dwarf star forming systems with GALEX and IRAS. 
Although significantly smaller than the SDSS, the dataset we have been building up over 
the last 10 years with data taken over a large stretch of the electromagnetic spectrum for 
few thousand galaxies in the local universe (worldwide available from the site GOLDMine; 
Gavazzi et al. 2003) turns out to be appropriate for the purposes of the present investigation. 
It includes drift-scan mode integrated spectra, narrow band Ha and broad band optical and 
near-infrared imaging for a volume limited sample of nearby galaxies in and outside rich 
clusters. The combination of GALEX and IRAS observations with these ancillary data allows 
us to study the dust attenuation properties in a sizable sample of normal star forming galaxies 
not suffering from the aperture bias and to compare observations with model predictions. 

In this first paper we investigate the relations between dust attenuation and global 
galaxy properties and compare them with the ones observed in starburst galaxies. The aim 
of this work is to provide some empirical relations based on observable quantities (thus model 
independent) suitable for deriving dust attenuation corrections when far infrared data are 
not available. For this reason all relations obtained throughout this paper will be given 
as a function of Ltir/ Lpuv, the observable that we consider the best dust attenuation 
indicator. We choose not to transform Ltir/ Lpuv into a (model dependent) estimate of the 
far ultraviolet extinction A{FUV), leaving the reader free to choose his/her preferred dust 
model (i.e. Meurer et al. 1999; Buat et al. 1999, 2002, 2005; Gordon et al. 2000; Panuzzo 
et al. 2003; Burgarella et al. 2005, Inoue et al. in preparation). Throughout the paper we 
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assume that quantities are related linearly and residual plots are presented in order to test 
the validity of this hypothesis. Moreover, since wc arc looking for new recipies to estimate 
the LtirI Lpuv ratio, this quantity has to be considered as the dependent variable, implying 
the use of an unweighted simple linear fit to estimate the best fitting parameters (Isobe et al. 
1990). 

A forthcoming paper will be focused on the comparison between models and observations 
hoping to gain a better understanding on the physics of dust attenuation and to study new 
recipes useful to convert Ltir/ Lpuv into A{FUV). 

2. The Data 

2.1. The optically-selected sample 

The analysis presented in this work is based on an optically selected sample of late-type 
galaxies (later than SOa) including giant and dwarf systems extracted from the Virgo Cluster 
Catalogue (VCC, Binggeh et al. 1985) and from the CGCG catalogue (Zwicky et al. 1961). 
The data include ~ 300 square degrees covering most of the Virgo, Abelll367 and Abell262 
clusters, the southwest part of the Coma cluster and part of the Coma-Al367 supercluster 
(llh30m < R.A. < 13h30m; 18° < decl. < 32°) observed in spring 2004 as part of the All-sky 
Imaging Survey (AIS) and of the Nearby Galaxy Survey (NGS) carried out by GALEX in 
two UV bands: FUV (Aeff = 1530A, AA = 400A) and NUV (Aeg = 2310A, AA = lOOOA). 
Details of the GALEX instrument and characteristics can be found in Martin et al. (2005) 
and Morrissey et al. (2005). Our sample has the quahty of being selected with the criterion 
of optical completeness. All galaxies brighter than a threshold magnitude are selected in 
all areas. In Coma-Al367 supercluster and A262 cluster all galaxies brighter than nip— 1^.1 
were selected from the CGCG catalogue (Zwicky et al. 1961). The Virgo region contains 
all galaxies brighter than mp= 18 from the VCC catalogue (Binggeli et al. 1985). We thus 
consider our sample an optically selected, volume limited sample. 

We include in our analysis all late-type galaxies, detected in both NUV and FUV GALEX 
bands and in both 60 /xm and 100 IRAS bands (157 objects). Whenever available, we 
extracted UV fluxes from the deep NGS images, obtained with a mean integration time of ~ 
1500 sec, complete to niAB ~ 21.5 in the NUV and FUV. Elsewhere UV fluxes have been ex- 
tracted from the shallower AIS images 70 sq. degrees), obtained with a mean integration 
time of ~ 100 sec, complete to itiab ~ 20 in both the FUV and NUV bands. All UV images 
come from the Internal Data Release vl (IRl.O). UV fluxes were obtained by integrating 
GALEX images within elliptical annuli of increasing diameter up to the optical B band 25 
mag arcsec"^ isophotal radii, consistently with the optical and near-lR images. Independent 
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measurements of the same galaxies obtained in different exposures give consistent pfiotomet- 
ric results within 10 % in the NUV and 15% in the FUV in the AIS, and a factor of ~ two 
better for bright (NUV <16) galaxies. The uncertainty in the UV photometry is on average 
a factor of ~ 2 better in the NGS than in the AIS, particularly for faint objects. The typical 
uncertainty in the IRAS data is 15% (Boselli et al. 2003). 

UV and far-infrared data have been combined to multifrequency data. These are optical 
and near-IR H imaging (mostly from Gavazzi et al. 2000, 2005; BoseUi et al. 2003), optical 
drift-scan spectra (Gavazzi et al. 2004; Gavazzi et al. in prep.) and Ho; imaging (Boselli & 
Gavazzi 2002; BoseUi et al. 2002a; Gavazzi et al. 1998, 2002b; Iglesias-Paramo et al. 2002; 
Gavazzi et al. in prep.), great part of which are available from the GOLDMine galaxy 
database (Gavazzi et al. 2003) (http:/ /goldmine. mib.infn.it). From the 157 galaxies selected 
we exclude Active Galactic Nuclei (AGN). AGNs have been selected using either the classi- 
fication provided by NED, if available, or by inspection to the integrated spectra of Gavazzi 
et al. (2004): we exclude galaxies with \og{[Oin]/Hp) > 0.61/(log([A^JJ]/i/a) - 0.05) + 1.3 
(Kauffmann et al. 2003). This criterion reduces the sample to 128 galaxies, spanning a range 
of six magnitudes in B band (-22< Mb <-16) and of three orders of magnitude in mass^ 
(9 < M < 12 Mq). Unfortunately ancillary data are not available for all galaxies observed 
by GALEX, we thus further divided the data in two subsamples. Sixty six galaxies in the 
primary sample have all the necessary complementary data (e.g. Ha photometry, Ha/H/? 
ratio, metallicity, H-band photometry; see Gavazzi et al. 2000, 2002a,b, 2004 for the selection 
criteria adopted in each survey). The remaining 62 galaxies form the secondary sample. We 
cannot exclude a possible contamination of AGN in the secondary sample, since no spectra 
are available for these objects. In all figures objects belonging to the primary sample will be 
indicated with filled circles while the secondary sample as empty circles. Since only galaxies 
belonging to the primary sample are present in all the plots analyzed in the presented work, 
all correlations will be quantified using only the primary sample. Data from UV to near-IR 
have been corrected for Galactic extinction according to Burstein Sz Heiles (1982). 
We assume a distance of 17 Mpc for the members of Virgo Cluster A, 22 Mpc for Virgo 
Cluster B, and 32 Mpc for objects in the M and W clouds (Gavazzi et al. 1999). Members of 
the Cancer, A1367, and Coma clusters are assumed to lie at distances of 65.2, 91.3, and 96 
Mpc, respectively. Isolated galaxies in the Coma supercluster are assumed at their redshift 
distance, adopting Ho = 75 km s^^ Mpc^^. 



Computed using the relation between Lh and M by Gavazzi et al. (1996) 
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2.2. The steirburst sample 

In order to compare the properties of our sample with starbursts, we compile a dataset 
of starburst galaxies observed by lUE from the sample of Calzetti et al. (1994). We consider 
29 galaxies, excluding AGNs and galaxies that have not been observed by IRAS at 60 or 100 
/im. Complementary data such as FIR, Ha fluxes and Balmer decrements are taken from 
Calzetti et al. (1995), metaUicities come from Heckman et al. (1998) and H-band photometry 
(available only for 18 galaxies) from (Calzetti 1997). Excluding the far infrared fluxes, all 
these quantities are obtained within an apertures of ~ 20 x lOarcsec^, consistent with lUE 
observations Calzetti et al. (1994). Thus we stress that aperture effects could strongly affect 
any comparison with normal galaxies for which all data are homogeneously integrated values. 
First of all, if the UV emission is more extended than lUE field of view the Ltir/ L},-uv ratio 
is overestimated^. In addition, even when physical quantities are obtained in the same lUE 
apertures, the presence of age and metallicity gradients in galaxies makes not trivial any 
comparison with the integrated values obtained for normal star forming galaxies (Kewley 
et al. 2005). All the observables, but the ultraviolet spectra slope /3, are cahbrated in a 
consistent way with our sample of normal galaxy. The ultraviolet spectral slope of starbursts 
is obtained by fitting lUE spectra (Calzetti et al. 1994), while for GALEX observations it 
comes from the FUV-NUV color index (see next Section). However, as shown by Kong et al. 
(2004), these two calibrations are consistent each other and do not introduce any systematic 
difference between the two samples. 



3. The Ltir/Lfuv — P relation for normal stair-forming galztxies 

Meurer et al. (1999) have shown that the ratio of far infrared to far ultraviolet luminos- 
ity tightly correlates with the UV colors of starburst galaxies. This relation, known as the 
infrared excess-ultraviolet (IRX-UV) relation, is often presented as f3 vs. Ltir/ Lpjjy rela- 
tion. As discussed in the introduction, we will refer throughout this paper to the Ltir/Lfuv 
ratio as the best indicator of UV dust attenuation and we will calibrate on it all the follow- 
ing relations. In order to determine the dust emission, we compute the total infrared flux 
emitted in the range 1-1000 jim, following Dale et al. (2001): 

log(/T7i?) = log(/F7fl) + 0.2738 - 0.0282 x log(^) + 

JlOO 



^However Meurer et al. (1999) argued that the majority of UV flux for their starburst sample Hes within 
the lUE aperture 
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+0.7281 X log(^)2 + 0.6208 x log(^)=' + 

/lOO /lOO 

+0.9118 X log(^)^ (1) 

JlOO 

where Jfir is the far-infrared flux, defined as the flux between 42 and 122 /im (Helou et aL 
1988): 

fpiR = 1.26 X (2.58 X /eo + /loo) x IQ-^^ [Wm'^] (2) 

and /eo and /loo are the IRAS fluxes measured at 60 and 100 /im (in Jansky). The total 
infrared luminosity is thus: 

Ltir = ^TvD^fTIR (3) 

The (3 parameter as determined from GALEX colors is very sensitive to the galaxy star for- 
mation history (see for example Calzetti et al. 2005). For this reason we assume throughout 
this paper (3 as defined by Kong et al. (2004): 

^~ -0.182 
= 2.201 X {FUV - NUV) - 1.804 (4) 

where fpuv and Jnuv are the near and far ultraviolet observed fluxes respectively (in 
erg cm s -1 A-i), and FUV and NUV are the observed magnitudes. 

The relationship between the ratio of total infrared luminosity (Ltir) obtained from (1) 
to the far-ultraviolet fluxes and the UV spectral slope (3 (or the FUV-NUV color) for our 
sample of nearby star forming galaxies is given in Fig.l. Several functional forms of the 
Ltir/ Lpuv — P relation can be found in the literature (i.e. Meurer et al. 1999; Kong et al. 
2004); we simply adopt a linear fit: log^LTm/ Lpuv) = a x (3 + b. This functional form is 
consistent with other previously proposed for /3 > — 2, while it diverges for /3 < — 2. Since 
the majority of normal and starbursts galaxies have (3 > —2 our choice is justified. This 
represents the simplest and less parameter dependent way to study the relation between two 
quantities.^ We find a strong correlation (Spearman correlation coefficient Vg ~0.76 for the 
primary sample and Vg ~0.65 for the secondary sample, both corresponding to a probabil- 
ity P{rs) >99.9% that the two variables are correlated) between the total infrared to far 
ultraviolet ratio and the spectral slope, but significantly different from the one observed for 
starburst galaxies (dashed line in Fig.l; Meurer et al. 1999). A test rejects at a confidence 
level higher than 99.9%, that the two samples follow the same relation. The best linear fit 



•^We tested this hypothesis fitting our data with functional forms similar to the ones proposed by Meurer 
et al. (1999) and Kong et al. (2004): no significative improvement in the scatter of this relation is obtained. 
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for our primary sample (solid line in Fig.l) is: 

log(^^) = (0.70 ± 0.06) xp+ (1.30 ± 0.06) (5) 

The uncertainty in the estimate of the Ltir/ Lpuv using equation (5) is ~ 0.26 ± 0.02 dex 
for the primary sample but it increases to ~ 0.35 ±0.03 dex, if we consider the whole sample 
(e.g. primary and secondary samples), consistent with the mean uncertainty observed for 
starburst galaxies (Meurer et al. 1999). A large contribution (~ 0.21 ± 0.02 dex) to the 
observed scatter in Eq.(5) is due to the uncertainty on the estimate of Ltir/Lfuv and /3. 
This result confirms once more that the Ltir/ Lpuy — {3 relation for normal galaxies deviates 
from the one observed for starbursts, as pointed out by previous studies of nearby galaxies 
(i.e. BeU 2002; Kong et al. 2004; Boissier et al. 2005; Buat et al. 2005; Seibert et al. 2005; 
Burgarella et al. 2005, Boissier et al. in prep.) and individual HII regions in nearby galaxies 
(Calzetti et al. 2005). 



3.1. The dependence on the birthrate parameter 

What physical mechanisms drive the difference observed in the Ltir/ Lpuv ~ P between 
normal star forming galaxies and starbursts? Recently Kong et al. (2004) interpreted the 
offset as an effect of the different star formation history experienced by galaxies and proposed 
that the distance from the starburst IRX-UV can be predicted using the birthrate parameter 
h (e.g. the ratio of the current to the mean past star formation activity, Kennicutt et al. 
1994). In order to test if the perpendicular distance dg from the Ltir/Lfuv — (3 relation 
for starbursts correlates with the star formation history of normal galaxies, we compute the 
birthrate parameter following BoseUi et al. (2001): 

SFRtoil - R) 



LH(Mt,t/LH){l-DM,,,,t) 

where R is the fraction of gas that stellar winds re-injected into the interstellar medium 
during their hfetime (~ 0.3, Kennicutt et al. 1994), to is the age of the galaxy (that we 
assume ~12 Gyr), DMcont is the dark matter contribution to the Mtot/Ln ratio at the 
optical radius (assumed to be 0.5; BoseUi et al. 2001). We compute the H-band luminosity 
following Gavazzi et al. (2002a): 

log Lh = 11.36 -QAx H + 2x log(L') [L©] 



where D is the distance to the source (in Mpc), and the SFR from the Ha luminosity 
(corrected for [Nil] contamination and for dust extinction using the Balmer decrement, see 
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Appendix A) following Boselli et al. (2001): 

Fig. 2 shows the relation between the birthrate parameter (eq.6) and the distance from the 
Ltir/ Lpuv — P relation for starburst galaxies. The two quantities are correlated (r^ ~0.40, 
corresponding to a correlation probability P{rs) ~99.8%) but with a large scatter. Given 
the value of observational uncertainties, it is not worth trying to use the observed trend 
to reduce the dispersion in the Ltir/ Lpjjy — {3 relation for normal galaxies. This result 
confirms that part of the dispersion in the Ltir/ Lpuv ~ P relation for normal star forming 
galaxies appears an effect of the different star formation history experienced by galaxies, as 
proposed by Kong et al. (2004). 



4. The P-A{Ha) relation 

Calzetti et al. (1994) found a strong relationship between the ultraviolet spectral slope 
(5 and the Balmer decrement Ha/H/3. For our starburst sample these two quantities are 
correlated (r^ ~0.81) as follows (see also blue stars in Fig. 3): 

p = (0.75 ± 0.10) X A{Ha) - (1.80 ± 0.13) (8) 

This empirical relation was used by Calzetti et al. (1994) to deduce an attenuation law (the 
Calzetti law), often applied to high redshift galaxies (i.e. Steidel et al. 1999; Glazcbrook 
et al. 1999). Contrary to the Ltir/ Lpuv — 13 relation the Calzetti law has not yet been 
tested for a sample of normal star forming galaxies. Buat et al. (2002) showed that for 
normal star forming galaxies the attenuation derived from the Calzetti law is ~0.6 larger 
than the one computed from FIR/UV ratio. This result has been recently confirmed by 
Laird et al. (2005) on star forming galaxies at 2; ~1. In order to check the Calzetti law on 
our sample we use the measure of the Hq;/H/3 described in Appendix A. Fig. 3 shows the 
relation between /3 and A{Ha) for our sample (empty and filled circles). For the primary 
sample we obtain ~0.58 (-P(rs) >99.9%) and: 

/5 = (0.37 ± 0.07) X A{Ha) - (1.15 ± 0.08) (9) 

flatter than for starburst galaxies (sec Fig. 3). At low A{Ha) normal galaxies show on average 
a less steep ultraviolet spectral slope than starbursts. In addition normal galaxies with the 
same value of P span a range of ~ 1 mag in A{Ha). At higher attenuation the two samples 
appear consistent. Our result suggest that the Calzetti law cannot be applied to normal 
galaxies. On the contrary, the relation between (3 and A{Ha) for normal galaxies, could be 
used to obtain a new attenuation law. 
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5. Relations between dust attenuation indicators and global properties. 

5.1. Metallicity 

Heckman et al. (1998) have shown that the ultraviolet spectral slope and metaUicity of 
starbursts are well correlated. To determine the metal content of our galaxies we average five 
different empirical determinations based on the following hne ratios: R23 = ([OII]A3727 + 

[OIII]A4959,5007)/H/3 (Zaritsky et al. 1994; McGaugh 1991), [NII]A6583/[01I]A3727 (Kew- 
ley k Dopita 2002), [NII]A6583/Ha (van Zee et al. 1998) and [OIII]A5007/[N11]A6583 (Dutil 
& Roy 1999). The mean uncertainty in the abundances is O.lOdex. In Fig. 4 we study the 
relationship between the gas metallicities and the LtirJ Lpuv ratio (left-panel) and j3 (right- 
panel) for normal star forming and starburst galaxies. For normal galaxies the Lxm/ Lpuv 
ratio correlates (r^ ~0.59, P{rs) >99.9%) with the gas abundance: 

log(^^) = (1.37 ± 0.24) X 12 + \og{0/H) - (11.36 ± 2.11) (10) 

with a dispersion of ~ 0.35 ± 0.03 in Xog^LTiRj Lpuy). As for the Ltir/ Lpuv — P relation 
normal galaxies differ from starbursts. At comparable metallicity normal galaxies show a 
lower Ltir/ Lpuv (lower attenuation) than starbursts, in agreement with the recent result by 
Boissier et al. (2004) who studied radial extinction profiles of nearby late-type galaxies using 
FOCA and IRAS observations. Unexpectedly we find however that normal star forming 
galaxies follow exactly the same (significant, ~ 0.58, P{rs) >99.9%) relationship between 
metallicity and ultraviolet spectral slope (5 determined for starbursts by Heckman et al. 
(1998) (see right panel of Fig.4). This might indicate that even though a normal and a 
starburst galaxy with similar gas metallicity have similar UV spectral slopes, they suffer 
from a significantly different dust attenuation, perhaps suggesting a different dust geometry 
(Witt & Gordon 2000). However we stress that this effect might occur due to aperture 
effects in the lUE data: while {3 is not significantly contaminated by aperture effects, the 
Ltir/ Lpuv ratio could be overestimated producing the observed trend (the total infrared 
luminosity is obtained by integrating the IRAS counts over the full galaxy extension, while 
the ultraviolet one is taken from lUE's significantly smaller aperture 20 x 10 arcsec^). This 
idea could be supported by the correlation (vg ~ 0.49, P{rs) >99.9% see Fig.5) observed 
between the starbursts' optical diameters and the Ltir/ Lpuv ratio, completely absent in 
our sample of normal galaxies (r^ ~ 0.006, P{rs) ~25%). GALEX observations of starburst 
galaxies will rapidly solve this riddle. 

We also checked the dependence of the dust attenuation on gas to dust ratio. Unfortunately, 
given the large errors on the estimate of the dust mass, we only obtain a weak relationship 
between the two quantities (see Appendix B). 
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5.2. Luminosity 



Since it is well known that the metallicity of normal galaxies strongly correlates with 
galaxy luminosity (e.g. Skillman et al. 1989; Zaritsky et al. 1994) and mass (e.g. Tremonti 
et al. 2004), it is worth considering the correlation between attenuation and galaxy lumi- 
nosity. Fig. 6 shows the relationships between the dust attenuation indicators Ltir/ Lpuy 
and /3 and the H-band luminosity. The infrared to far ultraviolet ratio correlates (r^ ~0.49, 
P{rs) >99.9%) with the total H-band luminosity: 



The dispersion of this relation is ~ 0.39 ± 0.03 in \og{LTiR/ Lpuv)- Since the H-band lumi- 
nosity is proportional to the dynamical mass (Gavazzi et al. 1996), this implies a relationship 
between dust attenuation and dynamical mass. Also in starbursts the total H-band lumi- 
nosity is correlated (r^ ~0.37, Pij-g) ~99.5%) with the Ltir/ Lpuv ratio and the great part 
of starbursts appear offset (to 99% confidence level) from the relation of normal galaxies. 
On the contrary, no difference is observed between the two samples in the i3-Lh plot, in 
agreement with what observed for metallicity. Finally Fig. 7 shows the relation between the 
bolometric luminosity (Ltir + Lpuv) and the dust attenuation, computed assuming that the 
UV emission is absorbed by dust and emitted in the far infrared. The correlation coeffi- 
cient {vg ~0.31, -P(rs) ~98%) indicates that the two quantities correlate, as for starburst 
galaxies (Heckman et al. 1998). This is not the case if we examine the relation between the 
ultraviolet spectral slope f3 and the bolometric luminosity (Fig. 7 right panel): while there is 
no correlation (r^ ~ 0.002, P{rs) ~20%) for our sample of normal galaxies, a clear relation 
{Ts ~ 0.68, P{rs) >99.9%) holds for starbursts. Starbursts with higher bolometric luminos- 
ity (high TIR emission) show lower ultraviolet slope, consistent with the idea that high TIR 
emission corresponds to high attenuation (low /3). 



Wang & Heckman (1996) interpreted the increase of dust attenuation with rotational 
velocity (or mass) as due to the variations in both the metallicity and surface density of 
galactic disk with galactic size. Fig. 8 shows the variation of the effective H-band surface 
brightness (defined as the mean surface brightness within the radius that contains half of the 
total galaxy hght) and the dust attenuation. The two quantities are strongly anti-correlated 



log(^^) = (0.34 ± 0.10) X log(^) - (2.66 ± 0.88) 




(11) 



5.3. Surface brightness 



(r, ~-0.63, P{rs) >99.9%): 




'FUV 



) = (-0.28 ± 0.04) X iie{H) + (5.92 ± 0.81) 



(12) 
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with a scatter of ~ 0.34 ± 0.03 in log{LTiR/ Lpuv)' ~1.2(T lower than the value obtained 
for H-band luminosity and consistent with the one obtained for the gas metallicity. Un- 
fortunately in this case we cannot compare the behaviour of normal galaxies with the one 
of starbursts due to the lack of an estimate of fi^ for the starbursts. Docs this relation 
indicate that UV dust extinction depends on the thickness of stellar disk, or does it follows 
from the correlation between attenuation and star formation surface density? To attack this 
question we determine the SFR density (defined as the ratio between the SFR determined 
from HcK (eq.7) and optical galaxy area). Fig. 9 shows the relation between the SFR density 
and log{LTiR/LFuv)- The two quantities are correlated {vg ~0.44, P{rs) >99.9%) with a 
dispersion of ~ 0.39 ± 0.03 in \og{LTiR/ Lpuv), ~1-2(J larger than the one observed for the 
mean H-band surface brightness^. Since the contribution of observational uncertainties to 
the scatter in the two relations is ~ the same (0.18 ± 0.02), our result might suggest that 
the UV attenuation is primarily correlated with the thickness of stellar disk, supporting 
the hypothesis of Wang & Heckman (1996) that both gas metallicity and star surface den- 
sity are directly connected with the physical properties of dust (i.e. quantity and spatial 
distribution) . 



5.4. LHa/Lpuv ratio 

Buat et al. (2002) suggested that the Lua/ Lpuv ratio could be another potential atten- 
uation indicator but they found a scattered correlation between Lho/Lfuv and A{FUV), 
confirmed by Bell (2002). This correlation is expected since both Ha and UV emission are 
star formation indicators. The Ha luminosity comes from stars more massive than 10 Mq 
and it traces the SFR in the last < 10^ yr while the UV luminosity comes from stars of 
lower mass (M> 1.5 Mq) and it can be used as an indicator of the SFR in the last 10^ 
yr. This means that under the condition that the star formation is approximately constant 
in the last ~ 10^ yr the ratio Lho/Lfuv (corrected for attenuation) should be fixed. Thus 
the ratio between the extinction corrected L^^ and the observed Lpuy should be a potential 
attenuation indicator. In Fig. 10 we analyze the relationship between the dust attenuation 
and the Lna/Lpuv ratio, where Lna is the Ha luminosity corrected for dust attenuation 
using the Balmer decrement and for the contamination of [Nil] . The two quantities turn out 
to be strongly correlated (r^ ~0.76, P{rs) >99.9%): 

log(^^) = (0.84 ± 0.07) X log(^^) - (0.59 ± 0.12) (13) 

Lfuv Lpuv 



^The difference between the two relations does not change if instead of the half-light radius, we use the 
total radius to estimate /Ue(-ff) 
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The dispersion around this relation is ~ 0.24 ±0.02 in log{LTiR/LFuv), consistent with the 
one observed for the log^LTm/Lpuv) ~ P relation. The high correlation and low scatter be- 
tween the two quantities is expected since the two variables arc mutually related: the FUV 
luminosity appears in both axes and Ltir and are known to be correlated (Kewley et al. 
2002), explaining why in the left panel of Fig. 10 starbursts and normal galaxies show the 
same trend. The right-panel of Fig. 10 shows the relation between the ultraviolet slope and 
the Lna/ Lpuv ratio. In this case starbursts and normal galaxies behave differently: at any 
given (3 starbursts have an higher L Ha/ Lfuv than normal galaxies, consistent with what 
expected for galaxies experiencing a burst of star formation (Iglesias-Paramo et al. 2004). 
A secure determination of the Balmer decrement for large samples is still a hard task, espe- 
cially at high redshift, thus we look for a relation similar to Eq.(13) using the observed Ha 
luminosity (Lf^a)- '^^^ L^a/LFUV and \og{LTiR/ Lfuv) ratios arc yet correlated (sec Fig. 11) 
but the correlation coefficient is lower than the previous case (r^ ~0.49, P{rs) >99.9%). The 
best linear fit gives: 

r T obs 

log(7^) = (1.10 ± 0.17) X log(-P^) - (0.59 ± 0.21) (14) 
Lfuv Lfuv 

with a mean absolute deviation of ~ 0.34 ± 0.03 (~3.3(7 higher than for Eq.l4). 

6. A cookbook for determining Ltir/ Lfuv ratio in optically-selected galatxies 

In this paper we investigated the relations between dust attenuation, traced by the 
Ltir/ Lfuv ratio, and other global properties of normal star forming galaxies. Furthermore 
we compared the dust attenuation in normal and starbursts galaxies using multiwavelength 
datasets. The amount of dust attenuation is found to correlate with the UV colors, gas 
metallicity, mass and mean surface density but, generally speaking, differently for normal 
and starburst galaxies. Determine whether this difference is real or is due to aperture ef- 
fects requires the analysis of GALEX observations for a sample of starburst galaxies. The 
dispersion in the Ltir/ Lfuv — P relation correlates with the birthrate parameter b, sug- 
gesting that the observed scatter is, at least partly, due to differences in the star formation 
history. These results stress that estimating the UV dust attenuation, and consequently the 
star formation rate of normal galaxies (at high redshift in particular) is highly uncertain 
(>50%) when rest- frame far infrared observations are not available. Moreover the sample 
selection criteria could strongly affect its properties, as recently pointed out by Buat et al. 
(2005) and Burgarella et al. (2005). They studied the dust attenuation properties and star 
formation activity in a UV and in a FIR selected sample, showing that the former shows 
correlations with global galaxy properties, such as mass and bolometric luminosity, that the 
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FIR selected sample does not. Their results stress that the dust attenuation properties are 
very heterogeneous and that Ltir/ Lpuv cannot be derived in a robust manner when far 
infrared observations are not available. 

However the present investigation has shown that among optically-selected samples of normal 
galaxies with no nuclear activity a number of empirical relations exists, allowing to derive 
the Ltir/ Lpuv ratio (and its uncertainty). Once the attenuation at UV is determined it 
can be transformed to any other A, only knowing the shape of the attenuation law and dust 
geometry (i.e. Calzetti et al. 1994; Gavazzi et al. 2002a; Boselli et al. 2003). 
In Table 1 we list all the relations, their associated r.m.s., mean absolute deviation from the 
best fit (m.a.d.)^ and the Spearman correlation coefficient. 

Before we proceed describing our recipes, we have to investigate whether the scatter in these 
relations is physical or is only driven by observational uncertainties. In the latter case, in 
fact, our cookbook would not be very useful, since it would be valid only for observations with 
the same uncertainties than our datasets. For H-band luminosity, H-band surface brightness, 
L'^^/ Lfuv ratio and metallicity the contribution of observational uncertainties to the ob- 
served scatter varies from ~ 18% (r.m.s.~ 0.17±0.02) for Lh to ~40% (r.m.s.~ 0.21 ±0.02) 
for 12 -I- \og{0/H) and L°^l^/ Lpuv'- even accounting for the contribution of measurements 
errors, the relative difference in the scatter of these relations does not change. On the con- 
trary this confirms that the relation involving Lh is the one with the highest "physical" 
dispersion, while for the other three relations the scatter is similar. 

The situation is worse for the relations involving (3 and the Lfja/ Lpuv ratio: the contribu- 
tion of observational errors is ~70-76% (~ 0.21 ± 0.02). Thus it is impossible to determine 
which of these two relations has the lowest scatter and represents the best way to estimate 
dust attenuation without far infrared observations. We can conclude that observational er- 
rors could account for the difference scatter observed in the relations involving (5 and the 
Lffa/Lpuv ratio, but not for the difference observed in all the other relations. Our results 
can thus be used to suggest different ways to correct for UV dust attenuation, 
la) The Ltir/ Lpuv — /5 relation still represents one of the best way to quantify dust atten- 
uation. The uncertainty in the value of \og{LTiR/ Lpuy) is ~ 0.26 ± 0.03. lb) If the UV 
spectral slope (3 is unknown but we know Lna (corrected for attenuation) we can obtain 
the ultraviolet attenuation using equation (13), with a r.m.s. of 0.24 ± 0.02. This relation 
is valid under the assumption that the star formation rate is approximately constant in the 
last Ri 10^ yr. 



^The mean absolute deviation is less sensitive to the contribution of outliers than the standard deviation. 

For a Gaussian distribution the mean absolute deviation (m.a.d.) is ^ ■\/2/7r x (r.m.s.), while it is lower 
(higher) for a heavier (lighter) tailed distribution. As shown in Table 1 the values obtained for r.m.s. and 
m.a.d. are consistent with the ones expected for a Gaussian distribution 
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Ila) If we know L'§^^, but no estimate of A(Hq;) is available, we can use Eq.(14) (rms~ 
0.34 ±0.03). 

lib) If neither f3 nor Ha luminosity are available we are left with the relations with H-band 
surface brightness^ (r.m.s.~ 0.34 ± 0.03) and, in the worse case, 
III) H-band luminosity (rms~ 0.39 ± 0.03 ). 

Summarizing, these relations allow us to estimate the value of the Ltir/ Lpuv ratio with an 
average uncertainties of~0.32 dex. This value corresponds approximately to a{A{FUV)) ~0.5 
mag, assuming \og{LTiR/ Lpuv) — 1 (the mean value for our sample) and using the model 
of Buat et al. (2005). This is the lowest uncertainty on the estimate of the Ltir/ Lpuv 
ratio in absence of far infrared observations. However we caution the reader that this value 
holds only for an optically-selected sample and that samples selected according to different 
criteria, especially FIR-selected, could contain higher dispersions. 
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A. Estimate of A{Hq.) 

The attenuation in the Balmer lines can be deduced from the comparison of the observed 
ratio Lfja/ Lufi with the theoretical value of 2.86 obtained for the recombination case B, an 
electronic density rif. < 10"^ cm^^ and temperature ~ 10^ K. The variation of this value 
with density its negligible and with temperature is <5% (in the range between 5000 K and 
20000 K, Caplan & Deharveng 1986). The underlying absorption was deblended from the 



^Since we need Ha flux to estimate metallicity, Eq.(lO) cannot be used in this case. 
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H/3 emission line using a multiple component fitting procedure. To do this the emission 
line is measured and subtracted from the spectra. The resulting absorption line is also 
measured with respect to a reference continuum. These two measurements are used as 
first guess in a fitting algorithm which fits jointly the emission and absorption lines to the 
reference continuum. For objects whose HI3 is detected in emission but the deblending 
procedure is not applied (no absorption feature is evident) a mean additive correction for 
underlying absorption equal to -1.8 in flux and -1.4 Ain EW is used. These values correspond 
to the fraction of the (broader) absorption feature that hes under the emission hne. We 
adopt a dust screen geometry and the Milky Way extinction curve (e.g. Kennicutt 1983; 
Calzetti et al. 1994). Whereas varying the extinction curves has neghgible effects in the 
visible, the dust screen assumption seems to under-estimate the extinction by ~0.2 mag 
compared with the amount deduced from the measurements of the thermal radio continuum 
(Caplan & Deharveng 1986; Bell & Kennicutt 2001). We do not apply any correction for 
Ha underlying absorption (Chariot & Longhetti 2001). However, since all the objects have 
EW{)i{a -\- [Nil]) > 3A, the underestimate in the value of A(Ha) is negligible. In fact no 
change (at a 99% significance level) is observed comparing the best fits obtained in this work 
and the ones obtained adding to the Ha the same fixed underlying absorption used for H/3 
when the underlying is not detected. We assume that the errors on A{Ha) are mainly due 
to the uncertainty on the H/3 flux. These errors represent in fact the lower limits because we 
do not account for the uncertainty introduced by the fitting of the lines. They range from 
0.01 to 0.43 mag and are found strongly anti-correlated with EW(H/9) (see Gavazzi et al. 
2004). Adopting the definition of the Balmer decrement as in Gavazzi et al. (2004): 



Cl{Hp) 



logij^ X 
(133 



(Al) 



Since the A{Ha) attenuation is: 



A{Ha) = 1.086 




) 



(A2) 



epa - 1 2.86 Lufj 



Prom (Al) and (A2) we obtain: 



A{Ha) = 1.086 



1 



X 0.33 X C1{H(3) /n(10) 



(A3) 



e/3a — 



and assuming a galactic extinction law {e^a — 1.47) we derive: 



A{Ha) = 1.756 x C1{H13) 



(A4) 



A{Ha) — 0.85 mag is obtained on average, consistent with previous studies (e.g. Kennicutt 
1983, 1992; Thuan & Sauvage 1992; Kewley et al. 2002). Eleven galaxies have H/3 undetected 
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in emission but the underlying stellar absorption is clearly detected. For them we derive a 
3 X crlower hmit to the H/3 flux {/hp) using (Gavazzi et al. 2004): 

fHf3 <Sx rms(45oo-480o) X Ha{HWHM) (A5) 

assuming that Ha and Hf3 emission lines have similar HWHM (Half Width Half Maximum). 
As shown in Eq.(Al) a change in the theoretical value of the Lna/Lnis ratio would only 
produce a small (<5%) constant over (or under) estimate of the ionized gas attenuation, 
thus leaving unchanged the shape and dispersions of the observed relations, only affecting 
the values of the best fitting parameters. 



B. Dust to Gas ratio 

The correlation between attenuation and metallicity can be interpreted assuming that 
the ultraviolet radiation produced by star forming regions suffers a dust attenuation increas- 
ing with the dust to gas ratio, which correlates with metallicity. (e.g. Issa et al. 1990; Inoue 
2003). In order to check this hypothesis we compute the dust to gas ratio following Boselli 
et al. (2002b). In normal galaxies the dust mass is dominated by the cold dust emitting 
above ~200 jim.. The total dust mass can be estimated provided that the 100-1000 //m 
far-IR flux and the cold dust temperature are known. Fitting the SEDs of normal galaxies 
with a modified Planck law z/^ B^iT^), with P = 2 (Alton et al. 2000), the total dust mass 
can be determined from the relation (Devereux & Young 1990): 

Mdust = C>5AL''(e"/^<^-* - 1) Mo (Bl) 

where C depends on the grain opacity, Sxis the far-IR flux at a given wavelength (in Jy), 
D is the distance of the galaxy (in Mpc), T^ust is the dust temperature, and a depends on 
A. Only IRAS data at 60 and 100 iim are available for our sample and, given the strong 
contamination of the emission at 60 /im by very small grains, the 60 to 100 /im ratio does not 
provide a rehable measure of Tdust (Contursi et al. 2001). T^ust determined by Alton et al. 
(1998) consistently with Contursi et al. (2001), seems to be independent of the UV radiation 
field, of the metallicity or of the total luminosity (Boselli et al. 2002b). Therefore we will 
adopt the average value Tdust = 20.8 ± 3.2 K for all our galaxies introducing an uncertainty 
of ~50% on the estimate of Mdust (equation (Bl)). We then estimate the dust mass of the 
sample galaxies using (Bl) with C = 1.27 MoJy^^Mpc^^, consistent with Contursi et al. 
(2001), and a=144 K for Sx = 5*100,^111 (Devereux Sz Young 1990). The determination of the 
dust to gas ratio, in a way consistent with that obtained in the solar neighbourhood, requires 
the estimate of the gas and dust surface densities, thus of the spatial distribution of dust 
and gas over the discs. Unfortunately only integrated HI and H2 masses are available for 
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our spatially unresolved galaxies. It is however reasonable to assume that the cold dust and 
the molecular hydrogen are as extended as the optical disc (Alton et al. 1998; Boselli et al. 
2002b). To determine the mean HI surface density we adopt (Boselli et al. (2002b)): 

log Shi = 20.92(±0.17) - 0.65(±0.11) x {def{HI)) cm-^ 

where def(HI) is the galaxy HI deficiency. Thus the dust to gas ratio is obtained from 
the ratio of the dust surface density to the sum of molecular and neutral hydrogen surface 
densities. In Fig. 12 we compare the relation between the Ltir/ Lpuv ratio (left panel) and 
(3 (right panel) with the dust to gas ratio. The gas to dust ratio barely correlates with the 
Ltir/ Lpuv ratio (R~0.38). Contrary to metallicity, we do not find a significant correlation 
(R~0.11) with the ultraviolet spectral slope. This is probably due to the high uncertainty 
in our estimate of M^ust consequent to assuming the same temperature for all our galaxies 
{M^ust oc e"/^'*""*, thus small errors (~15%) on T^ust propagate onto ~50% errors on Mdust)- 
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Table 1: Linear realtions useful to estimate the Ltir/ Lpuv ratio {log^LxiR/ Lpuv) = a x 
X + h). 



X 


a 


h 


m.a.d.^ 


rms^ 






0.70 ±0.06 


1.30 ±0.06 


0.20 ±0.02 


0.26 ±0.02 


0.76 


12 + \og{0/H) 


1.37 ±0.24 


-11.36 ±2.11 


0.26 ±0.02 


0.35 ± 0.03 


0.59 


Lh/Lq 


0.34 ±0.10 


-2.66 ±0.88 


0.29 ±0.03 


0.39 ±0.03 


0.49 




-0.28 ± 0.04 


5.92 ±0.81 


0.25 ±0.02 


0.34 ±0.03 


-0.63 


Lna/ Lpuv 


0.84 ±0.07 


-0.59 ±0.12 


0.19 ±0.02 


0.24 ±0.02 


0.76 




1.10 ±0.17 


-0.59 ±0.21 


0.27 ±0.02 


0.34 ±0.03 


0.49 



"Mean absolute deviation from the best fit 
''Standard deviation from the best fit 
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2^ -2 2 

Fig. 1. — Ratio of the total infrared to far ultraviolet luminosity as a function of the ul- 
traviolet spectral slope (lower x-axis) and the FUV-NUV color (upper x-axis). Open circles 
indicates our secondary sample while filled circles represent the primary sample. The dashed 
line represents the best linear fit to starburst IRX-UV relation. The solid line indicates 
the best bisector linear fit for our primary sample. The stars indicate the sample of lUE 
starbursts. Mean error bars for the plotted data are shown in the lower right corner, in this 
and subsequent figures. The residuals from the best linear fit for normal galaxies are shown 
in the bottom panel. 
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Fig. 2. — Relation between the birthrate parameter computed from the Ha emission, and 
the distance from the Ltir/ Lpuv ~ P relation for starbursts. The solid line represents the 
best linear fit. 
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Fig. 3. — The relation between the ultraviolet spectral slope (3 and the Ha attenuation 
obtained from the Balmer decrement. Symbols are as in Fig.l. Solid line represents the best 
linear fit to our primary sample (equation 9) while the dashed line indicate the best-fit for 
starburst galaxies obtained by Calzetti et al. (1994) (equation 8). Arrows indicate galaxies 
for which the value of A(Hq;) is a lower limit of the real value (i.e. H/3 undetected). The 
residuals from the best linear fit for normal galaxies are shown in the bottom panel. 
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Fig. 4. — Relation between gas metallicity and the Ltir/ Lpuy ratio (left) or [3 (right). 
Symbols are as in Fig.l. The solid lines show the best linear fit for our primary sample. The 
residuals from the best linear fits for normal galaxies are shown in the upper panels. 
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Fig. 5. — Relation between the galaxy size and the Ltir/ Lpuv ratio for starburst (left panel) 
and normal galaxies (right panel). Symbols are as in Fig. 1. Mean values and uncertainties 
in bins of 0.30 \og{Diameter) are given. 
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Fig. 6. — Relation between the H-band luminosity and the Ltir/ Lpuy ratio (left) or (3 
(right). Symbols are as in Fig. 1. The solid line shows the best linear fit for our primary 
sample. The residuals from the best linear fit for normal galaxies are shown in the upper 
panel. 
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Fig. 7. — Relation between the TIR+FUV luminosity and the LtirJ Lpuv ratio (left) or [3 
(right). Symbols are as in Fig. 1. 
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Fig. 8. — Relation between the mean H-band surface brightness (/ig) and the Ltir/ Lpuv 
ratio (left) or (5 (right). Symbols are as in Fig. 1. The solid line shows the best linear fit for 
our primary sample. The residuals from the best linear fit for normal galaxies are shown in 
the upper panel. 
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Fig. 9. — Relation between the star formation rate density and the Ltir/ Lpuv ratio (left) 
or /3 (right). Symbols are as in Fig. 1. The solid line shows the best linear fit for our 
primary sample.The residuals from the best linear fit for normal galaxies are shown in the 
upper panel. 
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Fig. 10. — Relation between the Ha and far ultraviolet luminosity and the Ltir/ Lpjjv ratio 
(left) or [3 (right). Symbols are as in Fig. 1. Ha luminosity is corrected for dust attenuation 
using the Balmer decrement, while the FUV flux is uncorrected. The solid lines show the 
best linear fit for our primary sample. The residuals from the best linear fit for normal 
galaxies are shown in the upper panels. 
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Fig. 11. — Relation between the observed Ha and far ultraviolet luminosity and the 
LtirI Lpuv ratio (left) or f3 (right). Symbols are as in Fig. 1. Ha luminosity is the 
observed value not corrected for dust attenuation. The solid lines show the best linear fit 
for our primary sample. The residuals from the best linear fit for normal galaxies are shown 
in the upper panels. 
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Fig. 12. — Relation between the gas to dust ratio and the LtirI Lpuv ratio (left) or (3 (right). 
Symbols are as in Fig. 1. The solid line shows the best linear fit for our primary sample. 



